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Abstract 

An approach, recently proposed to calculate the nucleon generalized par- 
ton distributions (GPDs) in a constituent quark model (CQM) scenario, in 
which the constituent quarks are taken as complex systems, is used to obtain 



helicity-dependent GPDs. They are obtained from the wave functions of the 



non relativistic CQM of Isgur and Karl, convoluted with the helicity-dependent 
GPDs of the constituent quarks themselves. The latter are modelled by using 
the polarized structure functions of the constituent quark, the double distri- 
bution representation of GPDs, and a phenomenological constituent quark 
form factor. The present approach permits to access a kinematical range 
corresponding to both the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi and the 
Efremov-Radyushkin-Brodsky-Lepage regions, for small values of the momen- 
tum transfer and of the skewedness parameter. In this kinematical region, the 
present calculation represents a prerequisite for the evaluation of cross sections 
relevant to deeply virtual Compton scattering. In particular, we have calcu- 
lated the leading twist helicity-dependent GPD H and, from our expressions, 
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its general relations with the non relativistic definition of the axial form factor 
and with the leading twist polarized quark density are consistently recovered. 



Typeset using REVTgX 



I. INTRODUCTION 



Generalized parton distributions (GPDs) [1-3] parametrize the non-perturbative hadron 
structure in hard exclusive processes (for a recent review, see, e.g., [4]). Their measurement 
would represent a unique way to access several crucial features of the structure of the nucleon, 
such as the Angular Momentum Sum Rule of the proton [3,5], as well as information on 
the structure of the proton in position space [6]. Therefore, relevant experimental efforts 
to measure GPDs, by means of exclusive electron Deep Inelastic Scattering (DIS) off the 
proton, are likely to take place in the next few years [7]. One of the most promising processes 
is represented by deeply virtual Compton scattering (DVCS) [2,3,8]. 

At twist-two, four helicity-even GPDs occur. Two of them are usually called unpolarized 
GPDs and labeled H and E. The first of them gives, in DIS kinematics, the usual quark 
density. The latter turns out to be a spin-flip distribution, in the sense that it implies a 
change of the spin of the target. The other two are usually called polarized or helicity- 
dependent GPDs. They are labeled H and E. The first of them yields in DIS kinematics 
the polarized quark density, the second is again a spin-flip distribution. The four of them 
enter the cross sections for the DVCS process, although the two non-spin-flip ones, H and 
H, give the dominant contribution. Reasonable estimates are necessary for the planning of 
experiments. 

For these quantities, several theoretical predictions have been already produced by using 
different descriptions of hadron structure: bag models [9], soliton models [10,11], light- 
front [12,13] and Bethe Salpeter approaches [14], and phenomenological estimates based on 
parametrizations of parton distribution functions (PDFs) [15,16]. Besides, an impressive 
effort has been devoted to study the perturbative QCD evolution [17] of GPDs, and the 
GPDs at twist three accuracy [18]. Chiral-odd GPDs deserve a special mention: while they 
have been shown to be in principle experimentally accessible [19], no theoretical estimates 
are available for them. 

We are interested in constituent quark model (CQM) calculations. In this framework, 
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in order to compare model predictions with data taken in DIS experiments, one has to 
evolve, according to perturbative QCD, the leading twist component of the physical structure 
functions obtained at the low momentum scale associated with the model, from this so called 
"hadronic scale", yu| [20,21], to the momentum scale of the data. Such a procedure has 
proven successful in describing the gross features of the standard PDFs by using different 
CQMs (see, e.g., [22]). Similar expectations motivated the study of GPDs in Ref. [23], 
where a simple formalism has been proposed to calculate the quark contribution to the 
GPDs from any non relativistic model. As an illustration, results for the GPD H q , obtained 
in the non relativistic (NR) model of Isgur and Karl (IK) [24] have been evolved from /Mq 
up to DIS scales, to next to leading order (NLO) accuracy. In Ref. [25] the same quark 
contribution has been evaluated, at /Xq, using the overlap representation of the GPDs [26] 
in light-front dynamics, along the lines developed in [27]. The same approach has been also 
applied to calculate the valence quark contribution to the helicity-dependent GPDs [28]. 
In all the papers dealing with CQM listed so far [23,25,28], only the Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) region of GPDs can be evaluated. To our knowledge, 
no CQM estimate of the helicity-dependent GPDs has been performed in the Efremov- 
Radyushkin-Brodsky-Lepage (ERBL) region. This is precisely the argument of the present 
paper, where an approach recently proposed for the unpolarized GPDs in a CQM scenario 
[29] , providing us with predictions also in the ERBL region, will be extended to the helicity- 
dependent GPDs. In Ref. [29], the procedure of Ref. [23] has been extended and generalized 
since, as a matter of fact, the latter, when applied to the standard forward case, has been 
proven to reproduce only the gross features of PDFs [22]. In order to achieve a better 
agreement with data, the approach has to be improved. In a series of papers, it has been 
shown that unpolarized [30] and polarized [31] DIS data are consistent with a low energy 
scenario, dominated by complex constituent quarks inside the nucleon, defined through a 
scheme suggested by Altarelli, Cabibbo, Maiani and Petronzio (ACMP) [32], updated with 
modern phenomenological information. The same idea has been applied to demonstrate the 
evidence of complex objects inside the nucleon [33], analyzing intermediate energy data of 
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electron scattering off the proton. A similar scenario, as old as the CQM itself [34], has been 
extensively used by other groups, starting form the concept of "valon" [35]. 

In Ref. [29] , the description of the unpolarized forward case of Ref . [30] has been general- 
ized and applied to the calculation scheme of Ref. [23]. As a result, more realistic predictions 
for the GPDs have been obtained and, at the same time, the ERBL region, not accessible 
before, has been explored. In particular, the evaluation of the sea quark contribution has 
become possible, so that GPDs can be calculated, in principle, in their full range of defini- 
tion. Such an achievement will permit to estimate cross-sections that are relevant for actual 
GPDs measurements, providing us with an important tool for planning future experiments. 

In here, the discussion is extended to helicity-dependent GPDs. The description of 
the polarized forward case of Ref. [31] is generalized to calculate helicity-dependent GPDs, 
following the path that lead to the results of Ref. [29] starting from the calculation of the 
unpolarized forward structure function in Ref. [30]. As in Ref. [29], the proposed approach 
will be applied here in a NR framework, which allows one to evaluate the GPDs only for 
small values of the 4- momentum transfer, A 2 (corresponding to A 2 <C m 2 , where m is the 
constituent quark mass) and small values also for the skewedness parameter, £. The full 
kinematical range of definition of GPDs will be studied in a followup, introducing relativity 
in the scheme. 

The paper is structured as follows. In the second section, after the definition of the main 
quantities of interest, a formula for the helicity-dependent current quark GPD H q , in a CQM 
with point-like constituents, will be derived. For the same quantity, in the third section, in a 
scenario dominated by composite constituent quarks, a convolution formula in terms of the 
corresponding quantity of the constituent quark and of a polarized off-diagonal momentum 
distribution, will be derived in impulse approximation (IA). Then, the helicity-dependent 
constituent quark GPDs are built in the fourth section, according to ACMP philosophy and 
using the Double Distribution (DD's) representation [16,36,37] of the GPDs. In the fifth 
section, results obtained by using CQM wave functions of the IK model and our helicity- 
dependent constituent quark GPDs will be shown. Conclusions will be drawn in the last 
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section. 



II. CALCULATIONS OF HELICITY-DEPENDENT GPDS IN MODELS WITH 



Let us consider a hard exclusive process, where the absorption of a high-energy virtual 
photon by a quark in a hadron target is followed by the emission of a particle to be later 
detected; finally, the interacting quark is reabsorbed back into the recoiling hadron. If the 
emitted and detected particle is, for example, a real photon, the so called Deeply Virtual 
Compton Scattering process takes place [2,3,10]. We adopt here the formalism used in Ref. 
[38]. Let us think of a nucleon target, with initial (final) momentum and helicity P(P') 
and s(s'), respectively. The GPDs H q (x, £, A 2 ) and E q (x, £, A 2 ) are defined through the 
expression 



where A = P' — P is the 4-momentum transfer to the nucleon, ip q is the quark field and M 
is the nucleon mass. It is convenient to work in a system of coordinates where the photon 
4-momentum, = (qo,q), and P = (P + P')/2 are collinear along z. The £ variable in 
the arguments of the GPDs is the so called "skewedness" , parametrizing the asymmetry 
of the process. It is defined by the relation £ = — n - A/2, where n is a light-like 4- vector 
satisfying the condition n • P — 1. As explained in [3], GPDs describe the amplitude for 
finding a quark with momentum fraction x + £ (in the Infinite Momentum Frame) in a 
nucleon with momentum (1 +C,)P and replacing it back into the nucleon with a momentum 
transfer A. Besides, when the quark longitudinal momentum fraction x of the average 
nucleon momentum P is less than — £, GPDs describe antiquarks; when it is larger than £, 
they describe quarks; when it is between — £ and £, they describe qq pairs. The first and 



POINT-LIKE QUARKS 




A • n7 5 
2M 



U(P,s) , (1) 
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second x intervals are commonly referred to as the DGLAP region and the third as the 
ERBL region [38], following the pattern of evolution in the factorization scale. One should 
keep in mind that, besides the variables x, £ and A 2 explicitly shown, GPDs depend, as the 
standard PDFs, on the momentum scale Q 2 at which they are measured or calculated. To 
make the presentation easier looking, we omit the latter dependence unless when specifically 
needed. The values of £ which are possible for a given value of A 2 are: 

< £ < v/^AVv^M 2 - A 2 . (2) 

The known constraints of H q (x, £, A 2 ) are: 

i) the so called "forward" limit, P' = P, i.e., A 2 = £ = 0, yields the usual polarized 
PDFs 

H q (x,0,0) = Aq(x) ; (3) 

ii) the integration over x, yields the contribution of the quark of flavor q to the axial 
form factor (ff) of the target A, called hereafter G^(A 2 ): 

JdxH q (x^,A 2 )=G q A (A 2 ) ; (4) 

iii) the polynomiality property [38], according to which the x-integrals of x n H q and of 
x n E q are polynomials in £ of order n + 1. 

In [23,29], the IA expression for the unpolarized H q (x, £, A 2 ), suitable to perform CQM 
calculations, was obtained. In a similar fashion, we derive next a formula for the helicity- 
dependent GPD H q (x,£, A 2 ) from the definition (1). 

Using the light-cone spinor definitions as given in the appendix B of [4], and defining: 
k + = (k + A; 3 )/v^2, k± = (fci, k 2 ), for the light-cone helicity combination s's — \\ — ++ one 
obtains 

Fl + (x, £, A 2 ) = Jl^im q (x, £, A 2 ) - -j*L^E q (x, £, A 2 ) , (5) 
so that, for £ 2 < 1: 
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Fl + (x, £, A 2 ) = H q (x, e, A 2 ) - e {^(x, £, A 2 ) + ^(x, £, A 2 )) + 0(£ 4 ) , (6) 



i.e. 



Fl + (x,Z,A*) = H q (x,Z,A*) + 0(e)- (7) 

According to this equation, in order to obtain the GPD H q (x, £, A 2 ) for £ 2 <C 1 one has to 
evaluate starting from its definition, Eq (1). Using light-cone quantized quark fields in 
the l.h.s., whose creation and annihilation operators, W(k) and b(k), obey the commutation 
relation 

{b(k'), b\k)} = (27r) 3 2k + 5(k' + - k + )5 2 (k ± - k' ± ) , (8) 
and using light-cone states normalized as 

(P'\P) = (2tt) 3 2P + 5(P ,+ - P + )5 2 (P ± - P|) , (9) 
one obtains, for x > £ [38] 



2^-^1(2^)5 
where 

(tf(k')b(k)) = J2 stgn(X)(P' + \b{((x - 0P + , k' ± )b x ((x + £)^ + , k±)\P+) , (11) 

A 

and V is a volume factor. Since we want to obtain only the quark contribution to H q , the 
only one which can be evaluated in a CQM with three valence, point-like quarks, we are 
interested here in the x > £ region. 

Eq. (10) can be written (see [29] where all the steps of a similar derivation are explicitely 
given) : 

*M*. {, A') = ^ / 2 ^ 2 ^ (£-(* + 0) <»'(* + A)«*» • (12) 

In a NR framework, states and creation and annihilation operators have to be normalized 
according to 
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(P'\P) = (2tt) 3 5(P'+ - P + )5(P' ± - P ± ) (13) 

and 

{b(k' + , k ± ), b\k + , k' ± )} = (27i) 3 5(k' + - k + )5(k' ± - k ± ) (14) 
respectively. As a consequence, one has to consider that [39]: 

\P) - V2P^\P) , (15) 



b(k) -> V2k+b(k + ,k ± ) , (16) 
so that in Eq. (12), in terms of the new states and fields, one has to perform the substitution 
(b\k + A)b(k)) = £ sign(X) (P' + \b[((x - £)P + , k± + &±)bx((x + k ± )\P+) - 



2yjl - £ 2 P + V2k>+2k+ x 

xY, sign(\) (P'\b[(k+ + A + ,k ± + A ± )b x (k + ,k ± )\P) , (17) 



where the relation 2^/1 — £ 2 P + = \ / 2P + \ / 2P' + has been used. 
Now, inserting Eq. (17) in Eq. (12), one gets 



1 r d 2 k ± dk + (k + 
V 1 2 v / FF+(2tt) 3 [p 1 
x J2 sign(X) (P'\b\(k + + A+ k ± + A ± )b x (k + , k ± )\P) = 



1 fd 2 k ± dk+Jk + . \ 
x £ sign(\)(P'\b[(k+ + A+,k L + A L )b x (k + ,k 1 )\P) + 0(e) ■ (18) 



The constituent quarks with mass m are taken to be on shell, k = y k 2 + m 2 , thus 

d 2 k ± dk + -> (tjpj dk , (19) 
and from Eq. (18) and Eq. (7), one gets 
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H q (x^,A 2 ) = J dksl^-ix + O^J x 
1 k+ 



x 



]T s^n(A) (P'|6l(A; + + A+, k ± + A ± )6 A (A;+, k ± ) \P) 



(2^Vk . 

' k 2 (k + A) 2 



> + 



n,(fc, fc + A) + O 



m 2 ' w? 



+ o(£ 2 ) 



In the last step, the definition of a (non-diagonal) polarized NR momentum distribution, 
h q (k,k + A), has been used, together with the fact that a NR momentum distribution 

— * 

describes the probability of finding a constituent of momentum A; in a given system up to 
terms of order [401. 

So far, we have shown that, in a NR CQM, the GPD H q (x, £,A 2 ) can be calculated, 
for £ 2 <C 1, k 2 m 2 and (fc + A) 2 <^ m 2 (which means, in turn, A 2 m 2 ), through the 
following expression: 

H g (x, C A 2 ) = J dk5 (~+-(x + £)) " 9 A ^ + A) + O ^, ^ j ■ (22) 

The above equation, corresponding to Eq. (8) in [23] if polarization is taken into account, 
permits the calculation of H q (x, £, A 2 ) in any CQM, and it naturally verifies some of the 
properties of the GPDs. In fact, the polarized quark density, Aq(x), as obtained by analyz- 
ing, in IA, the DIS off the nucleon (see, e.g., [41]), assuming that the interacting quark is 
on-shell, is recovered in the forward limit, where A 2 = 0, £ = 0: 

Aq(x) = H q (x, 0, 0) = J dkn q (k) 5 (x - , (23) 

so that the constraint Eq. (3) is fulfilled. In the above equation, n q (k) is the NR polarized 
momentum distribution of the quarks in the nucleon and, integrating Eq. (22) over x, one 
obtains the contribution of the quark of flavor q to the IA definition of the axial ff, so that 
Eq. (4) is verified. 

The polynomiality condition is formally fulfilled by the GPD defined in Eq. (22), al- 
though the present accuracy of the model, explicitly written in the latter equation, does not 
allow to really check polynomiality, due to the effects of the Munich Symmetry [42], which 
excludes 0(£) contributions. 
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With respect to Eq. (22), a few caveats, already listed in Ref. [29] for the unpolarized 
case, are necessary. 

i) One should keep in mind that Eq. (22) is a NR result, holding for £ 2 <C 1, under the 
conditions ? « m 2 , A 2 < m 2 . If one wants to treat more general processes, the NR limit 
should be relaxed by taking into account relativistic corrections. In this way, at the same 
time, an expression to evaluate E q (x,£, A 2 ) could be obtained. 

ii) The Constituent Quarks are point-like objects. 

iii) If use is made of a CQM, containing only constituent quarks (and also antiquarks 
in the case of mesons), only the quark (and antiquark) contribution to the GPDs can be 
evaluated, i.e., only the region x > £ (and also x < — £ for mesons) can be explored. In 
order to introduce the study of the ERBL region (— £ < x < £), so that observables like 
cross-sections, spin asymmetries and so on can be calculated, the model has to be improved. 

iv) In actual calculations, the evaluation of Eq. (22) requires the choice of a reference 
frame. In the following, the Breit Frame will be chosen, where one has A 2 = —A 2 and, 
in the NR limit we are studying, one finds \/2P + — > M. It happens therefore that, in the 
argument of the 5 function in Eq. (22), the x variable for the valence quarks is not defined 
in its natural support, i.e. it can be larger than 1 and smaller than £. Although the support 
violation is small for the calculations that will be shown here, it has to be reported as a 
drawback of the approach. 

The issue iii) will be discussed in the next sections, by relaxing the condition ii) and 
allowing for a finite size and composite structure of the constituent quark. 

III. HELICITY- DEPENDENT GPDS IN A CONSTITUENT QUARK SCENARIO 

In the standard forward case, the procedure described in the previous section, has been 
proven to be able to reproduce the gross features of PDFs [22]. In order to achieve a better 
agreement with data, the approach has to be enriched. 

In previous papers, it has been shown that unpolarized [30] and polarized [31] DIS 
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data are consistent with a low energy scenario dominated by complex systems of point-like 
partons. The forward parton distributions are therefore given by a convolution between 
constituent quark momentum distributions and constituent quark structure functions. The 
latter quantities were obtained by using updated phenomenological information in a scenario 
firstly suggested by Altarelli, Cabibbo, Maiani and Petronzio (ACMP) already in the sev- 
enties [32]. Following the same idea, in Ref. [29] a model for the reaction mechanism of an 
off- forward process, such as DVCS, where GPDs could be measured, has been proposed. As 
a result, a convolution formula giving the proton H q GPD in terms of a constituent quark off- 
forward momentum distribution, H qo , and of a GPD of the constituent quark q Q itself, H qoq , 
has been derived. In here, a similar program will be developed for the helicity-dependent 
GPD H q . 

It is assumed that the hard scattering with the virtual photon takes place on a parton of 
a spin 1/2 target, such as the proton, made of complex constituents. The scenario is depicted 
in Fig. 1 for the special case of DVCS, in the handbag approximation. One parton (current 
quark) with momentum k, belonging to a given constituent of momentum p, interacts with 
the probe and it is afterwards reabsorbed, with momentum k + A, by the same constituent, 
without further re-scattering with the recoiling system of momentum P R . We suggest here 
an analysis of the process which is quite similar to the usual IA approach to DIS off nuclei 
[40], recently applied also to nuclear GPDs [43]. 

In the class of frames chosen in section 2, and in addition to the kinematical variables, 
x and £ one needs a few more to describe the process. In particular, x' and for the 
"internal" target, i.e., the constituent quark, have to be introduced. These quantities can 
be obtained defining the "+" components of the momentum k and k + A of the struck parton 
before and after the interaction, with respect to P + and p + = \(p + p') + : 

k+ = (x + S)P + = (x' + £')P + , (24) 
(k + A)+ = (x - 0P + = {x' - C)P + ■ (25) 

From the above expressions, £' and x' are immediately obtained as 
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^ 2p+ 

x' = (27) 
and, since £ = — A + /(2P + ), if 5 = p + /P + , one also has 

^iorW (28) 

In order to derive a convolution formula in IA for H q , the procedure used for standard 
DIS will be adopted [40]. First, in Eq. (20), two complete sets of states, corresponding to 
the interacting particle and to the recoiling system, are properly inserted to the left and 
right-hand sides of the quark operator: 

H q (x,£, A 2 ) = (P'S\ £ {{PkSWs'miP'AKp's'l} 

P' w S> wP ',s> 

7^3 ^ - (* + O) \ E &J >A (* + + A+, fc ± + A ± )M^ fc ± ) 

E {l^«>bs)}{<PA|<ps|} 1^5); (29) 

Pr,Sr,p,s 

second, since eventually use has to be made of NR wave functions, a state \p) has to be 
normalized in a NR manner: 

tf\p) = (2ic)*5tf-d, (30) 
so that in Eq. (29) one has to perform the substitution: 

\P>^\J^\P>- (31) 
Since, using IA in the intrinsic frame of the target one has, for the NR states: 

{(P R S R \(ps\}\PS) = (P R S R ,ps\PS)(27r) 3 5 3 (P - P R - p)5 s ,s R s , 
a convolution formula is readily obtained from Eq. (29): 



tf g (:r,£,A 2 )^ W dE I d P\ P T 1 a\ + P ^P + K ' E )jH gog W, A 2 ) + O (e) = (32) 
qo J J \ Po{P + A Jo ? v 7 

= J2[ dE [ dp[P qo (p,P + A, E) + 0(^/M 2 , A 2 /M 2 )] 

x |^ 0? (x',£',A 2 )+O(£ 2 ) . (33) 



13 



In the above equations, the kinetic energies of the residual system and of the recoiling target 
have been neglected, and P qo (p,p + A, E) is a spin dependent, off-diagonal spectral function 
for the constituent quark in the target: 



P qo (p,p + A,E) = J-^Y, sign(s) (P'M\(P - p)S R , (p + A)s)((P - fiS R ,ps\PM) x 

x 5(E - E* R ) . (34) 



Besides, the quantity 



H qoq {x'^\A 2 ) = J dk5^—-(x' + i') 



x WfvT^ sl9<X) (M^ + + ^ + M + A ± )b q , x (k + ,k ± )\ P ) = 



1 



X WfVh^ s ^(A)(p , |&;, A ( A;+ + A+ '^ + ^)^(A; + ,^)b> (35) 

is, according to Eq. (20), the helicity-dependent GPD of the constituent quark q up to 
terms of order 0(£ 2 ), and in the above equation use has been made of Eqs. (27) and (28). 

The delta function in Eq (34) defines E, the removal energy, in terms of E R , the excitation 
energy of the recoiling system. 

Concerning Eqs (33-35), two comments are in order. 

The first concerns the accuracy of the actual calculations which will be presented. In the 
following, NR wave functions will be used to evaluate Eq. (33) , so that the accuracy of the 
calculation is of order O (p 2 /M 2 ,A 2 /M 2 ). In fact, if use is made of NR wave functions to 
calculate the non-diagonal spectral function, the result holds for p 2 , (p + A) 2 « M 2 . The 
same constraint can be written p 2 , A 2 << M 2 . 

Second, from the study of forward DIS off nuclei, it is known that, in going from a 
covariant formalism to one where use can be made of the usual NR wave functions, one has 
to keep the correct normalization [40] . This procedure leads to the appearance of the "flux 



factor", represented in Eq. (32) by the expression yp + (p + A) + /p (p + A) (which gives, 
in forward DIS, the usual p + /po term). This factor gives one at order O [p 2 jM 2 , A 2 /M 2 ), 
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which is the accuracy of the present analysis, and it will not be included in the calculation 
described below. If one wants to be predictive in more general processes at higher momentum 
transfer, the accuracy of Eq. (33) is not good enough any more and the calculation, in order 
to be consistent, has to be performed taking into account relativistic corrections. 
Eq. (33) can be written in the form 

A'tt 



Hq(x,£, A 2 ) = £ fdE f dpP qo (p,p + A)jH qoq (x', A 2 ) = 

qo J J S 

= E jf y j dE j dpP qo (p,P + A> [z - |) H qoq 



"A A 2 ) . (36) 



go 

Taking into account that 



z z 



z~l; = z-Ml + $-Z] = z + Z-^(l + $ = z + Z- ] ]^, (37) 



Eq. (36) can also be written in the form: 

H q (x, C A 2 ) = £ f -H qo (z, A*)H qoq (*, A 2 

qo Jx Z \Z Z 



(38) 



where 



H qo (z, C A 2 ) = jdEj dpP qo (p, ir+£)s(z + t-£^ . 



(39) 



One should notice that Eqs. (38) and (39) or, which is the same, Eq. (33), fulfill the 
constraint i) — in) previously listed. 

The constraint i), i.e. the forward limit of GPDs, is verified. In fact, by taking the 
forward limit (A 2 — > 0,£ — > 0) of Eq. (38), one gets the expression which is usually found, 
for the polarized parton distribution Aq(x), in the IA analysis of polarized DIS off polarized 
three-body systems made of composite constituents (see, i.e., [44] for the 3 He target) 

rl dz / X\ 

Aq(x) = H q (x, 0, 0) = J2 / —9M A fe - ■ ( 4 °) 

qo Jx Z \ZJ 

In the latter equation, 

g qo (z) = H qo (z, 0,0) = jdEj dpP qo (p, E)5 (z - (41) 
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is the light-cone momentum distribution of the constituent go in the nucleon, Aq qo (x) = 
H qoq (x, 0, 0) is the distribution of the quark of flavour q in the constituent q and P qo (p, E), 
the A 2 — ► limit of Eq. (36), is the one body spin dependent spectral function, labeled P\\ 
in Ref. [45,44]. If the system were a three-body nucleus, the constituents would be nucleons 
and Eqs. (38)-(41) would coincide with the corresponding ones obtained in Ref. [44] for the 
polarized 3 He target. 

The constraint ii), i.e. the x— integral of the GPD H q , is also naturally fulfilled. In fact, 
by x— integrating Eq. (38), one easily obtains: 

/ dxH q (x,£,A 2 ) = £ fdx f -H qo (z,£,A 2 )H qoq (-, -, A 2 

J J J Z \^ Z Z 

= W dx'H qoq (x',e, A 2 ) / dzH qo (z,^,A 2 ) = 
go J J 

= E^ (A 2 )G*(A 2 )=q(A 2 ). (42) 

go 

In the equation above, G q p (A 2 ) is the contribution, of the quark of flavour q, to the proton 
axial ff; G^ (A 2 ) is the contribution, of the quark of flavour q, to the constituent quark q 
axial ff; G* (A 2 ) is a proton axial "pointlike ff", which would represent the contribution 
of the quark q to the axial ff of the proton if q were point-like. G^°(A 2 ) is given, in the 
present approximation, by 

G«'(A 2 ) = jdEj dpP qo (p,p + A,E) = J dzH qo (z,Z,A 2 ) . (43) 

Eventually the polynomiality, condition Hi), is formally fulfilled by Eq. (33), although 
one should always remember that it is a result of order 0(£ 2 ), so that high moments cannot 
be really checked. 

Summarizing, an IA convolution formula has been derived, valid for any helicity- 
dependent GPD H q of a spin 1/2 hadron target made of three spin 1/2 constituents, Eq. (33) 
(or, which is equivalent, Eq. (38)), at order O (^2, ;pi;£ 2 )> i n terms of a spin-dependent, 
non-diagonal spectral function, Eq. (34), and of the helicity-dependent GPD H qoq of the 
internal target. Convolution formulae are usually obtained in describing quantities, such 
as structure functions, which have a probabilistic interpretation. This is the case of Eq. 
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(40). However, GPDs represent amplitudes and not probability densities. Our convolution 
equation for GPDs, obtained as the result of a mathematical procedure, also used by sev- 
eral authors in the context of the IA in nuclear GPDs [46,47,43], has a different meaning. 
It involves wave functions evaluated for different values of the momentum of one of the 
particles (one body off diagonal matrix elements). Besides, in the derivation, the simple 
structure of the final result, Eq. (38), is due to the restrictions imposed in the kinematics 
(A 2 <C m 2 ,£ 2 <C 1). Had these approximations, which are beyond IA, not been made, 
Eq. (38) would be more complicated. For example, a contamination from the GPD E q 
would be found. Therefore, the formulae we are discussing here are not strict convolution 
formulae in the conventional probabilistic sense. Finally, if £ = 0, the formula defining H 
becomes a strict convolution (cf. Eq. (5)) and, as a matter of fact, in this case GPDs have 
a probabilistic interpretation, being the Fourier transforms of the impact factor dependent 
parton distributions [6]. Our formalism is therefore consistent with the fact that convolu- 
tions, in the probabilistic sense, arise once we limit the expressions to quantities having a 
probabilistic interpretation. 

Let us now specify the formalism to the proton case and to CQM calculations. If the 
.E 1 - dependence of i.e., the E-dependence of z (cf. Eq. (28)) is disregarded in Eq. (33), so 
that the spin dependent, off-diagonal momentum distribution 



which can also be cast in the form of Eq. (38), with the light-cone off-diagonal momentum 
distribution given by 




(44) 



is recovered, Eq.(33) can be written in the form 




(45) 




(46) 
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which coincides with Eq. (22) and it is to be evaluated in a given CQM, for q Q = uq or do, 
while H qoq (j, |, A 2 ) is the constituent quark helicity-dependent GPD, which is still to be 
discussed and will be modelled in the next section. 

IV. A MODEL FOR THE HELICITY-DEPENDENT GPD OF THE 

CONSTITUENT QUARK 

The main issue remaining is the definition of H qoq (j, | , A 2 ), the helicity-dependent con- 
stituent quark GPD, appearing in Eq. (38). 

As it has been done in Ref. [29] for the unpolarized case, this quantity will be modelled 
to respect that its forward limit, recovers the constituent quark parton distributions. We 
adopt the picture proposed in [30,31] that the constituent quark is a complex system of 
point-like partons described by the scheme suggested by Altarelli, Cabibbo, Maiani and 
Petronzio (ACMP) [32]. In this scheme the constituent quarks are composite objects whose 
structure functions are described by a set of functions 4> qoq (x) that specify the number of 
point-like partons of type q which are present in the constituent of type go, with fraction x 
of the total momentum. We call these functions, generically, the structure functions of the 
constituent quark. 

The functions describing the nucleon parton distributions are expressed in terms of the 
independent 4> qoq (x) and of the constituent density distributions (q = u ,d ) as, 

q(x, Q 2 ) = W 1 -q (z, Q 2 )<P qoq (-, Q 2 ) , (47) 

where q labels the various partons, i.e., valence quarks (u v ,d v ), sea quarks (u s ,d s ,s), sea 
antiquarks (u, d, s) and gluons g. 

The different structure functions of the constituent quarks are derived from three very 
natural assumptions [32]: i) the point-like partons are QCD degrees of freedom, i.e. quarks, 
antiquarks and gluons; ii) Regge behavior for x — > and duality ideas; iii) invariance under 
charge conjugation and isospin. 
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These considerations define in the case of the valence quarks the following structure 
function 



iXA+j) (l-x) A ~ 

r(|)r(A) ^ 



^(x,Q 2 ) = ^luv^T- g ■ (48) 



For the sea quarks the corresponding structure function becomes, 



b qqs {x,Q 2 ) = -{l-x) D -\ (49) 

Ob 



and, in the case of the gluons, it is taken 

<j )qg (x,Q 2 ) = -(l-x) B - 1 . (50) 

Ob 

The last input of the approach consists of the definition of the scale at which the con- 
stituent quark structure is defined. We choose the so called hadronic scale /Xq [21,22]. This 
last hypothesis fixes all the parameters of the approach (Eqs. (48) through (50)). The 
constants A, B, G and the ratio C/D are determined by the amount of momentum carried 
by the different partons, corresponding to a hadronic scale of /1q = 0.34 GeV 2 , according to 
the parameterization of [21]. C (or D) is fixed according to the value of F 2 at x = 0.01 [32], 
and its value is chosen again according to [21]. We stress that all these inputs are forced 
only by the updated phenomenology, through the 2 nd moments of PDFs. The values of the 
parameters obtained are listed in [30]. 

The other ingredients appearing in Eq. (47), i.e., the density distributions for each 
constituent quark, represent the unpolarized equivalent of Eq. (22). 

In Ref. [31], the spin dependent structure function of the constituent quark have been 
built generalizing the analysis of Ref. [30], leading to Eqs. (48) - (50). In here, we review 
the procedure addressed in Ref. [31]. Let 

= q+(x,f4) - q-(x,n 2 ) (51) 

where ± label the quark spin projections and q represents any flavor. The AC MP approach, 
generalized to the polarized case, implies 

19 



(52) 

go ^ z j * 
where i = ± labels the partonic spin projections, j = ± the constituent quark spins and 
g = Uo,do. Using spin symmetry we arrive at 



Aq(x) = [ 

Jx 



1 dz 

x Z 



Au (z)A$ Uoq (-) + Ad (z)A$ doq (-)) 



(53) 



where Ag = <?o+ — qo- , an d 



^■^qaq ^qo+q+ ^qo+q- 



(54) 



Note at this point that the unpolarized case can be described in this generalized formalism 

as 



q{x) = q + (x) + q-(x) = ^ / — q (z)<f> qoq (-) , 

qo Jx Z Z 



(55) 



where 



®qoq — ®qo+g+ + 3*90+9- > 



(56) 



The description, reformulated in terms of the conventional valence and sea quark sepa- 
ration, yields 



Aq(x) = Aq v (x) + Aq s (x) = 



qo Jx Z 



, X 



X , 



A^(-) + A^ s (-) 



(57) 



Since AQ qoqv = A$ qoqv 5 qoq , i.e. any constituent has a leading parton with the same quantum 
numbers, one has 



Aq v ( x ) = f— Aq {z)A$ qoqv (-), 

Jx Z Z 



(58) 



qo Jx z z 



(59) 
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Introducing SU (6) (spin-isospin) symmetry as a simplifying assumption, gives 

A$ U()M = A$ dod (60) 

and 

A$ uod = A$ doU . (61) 

Furthermore they imply 

A$ U0Uv + A$ U0Us = A$ dQdv + A$ dods (62) 

and 

A$ Uods = A$ doUs . (63) 
By adding to these the following relations 

A$ U0Us = A$ doUs , (64) 

A$ d()ds = A$ Uods , (65) 
which are beyond SU (6) symmetry but quite reasonable, one obtains, 

A$ U0Us = A^ Us = A$ U()(is = A$ dods = A$ qoqs (66) 

and 

A$ UoUv = A$ dodv = A$ qoqv . (67) 

In this way the structure functions for the valence and for the sea are reduced to just 
two independent constituent structure functions and Eq. (59) simplifies to 

Aq s (x) = C -(Auo(z) + Ad (z))A<S> qoqs (-). (68) 

Jx Z Z 

The same argumentation applied to gluons implies 
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Ag(x) = / -(Au (z) + Ad (z))A% og (-) (69) 

Jx Z Z 

Thus the AC MP procedure can be extended to the polarized case by introducing three 
additional structure functions for the constituent quarks: A$ qoqv , A$ qoqs and A§ qog . 

In order to determine the polarized constituent structure functions we add some assump- 
tions which will determine them and those of the unpolarized case using a small number 
of parameters. They are: iv) factorization assumption: A$ cannot depend upon the quark 
model used, i.e, cannot depend upon the particular Ag ; v) positivity assumption: the pos- 
itivity constraint A$ < $ is saturated for x = 1. These additional assumptions determine 
completely the parameters of the polarized constituent structure functions. In fact, the 
QCD partonic picture, Regge behavior and duality imply that 

A% )f = AC f x a f(l-x) A f~ 1 (70) 

and < df < |, for all / = q v ,q s ,g, as defined by dominant exchange of the A 1 meson 
trajectory [48]. 

The positivity restriction, A$ < $, is a natural consequence of the probability inter- 
pretation of the parton distributions. The assumption that the inequality is saturated for 
x — 1, in the spirit of ref. [49], implies that AC/ = C/, the latter being the parameter fixed 
in the analysis of the unpolarized case, and therefore when x ~ 1 the partons which carry 
all of the momentum also carry all of the polarization, i.e., = 0. As is discussed in Ref. 
[31], this scheme has been changed only for the sea structure function, in order to obtain a 
reasonable description of the polarized sea. This led to the definition of a parameter, AC Qs , 
different from the one, C qs defined in the unpolarized case. From the point of view of the 
number of parameters, this is a minimal assumption, since it reduces the parameters of the 
polarized case to those of the unpolarized one, with the only exception of C qg . Lastly the 
exponents Af are also taken from the unpolarized case. To summarize the parameterization, 
let us stress that the main change between the polarized functions and the unpolarized ones 
comes only from Regge behavior and the only parameter which has been changed is the one 
related to the polarized sea distribution. 
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With all the above hypothesis the constituent quark functions become 

fi) = T^T^f' i 1 -^ 1 ( 71 ) 



A% qs (x,v 2 ) = AC s x a f(l-x) D - 1 (72) 



A$ 99 (x,/i 2 ) = Gx a f(l- x) B - 1 (73) 

where A, D, B, G are the parameters of the unpolarized case, while AC S is taken to reproduce 
the polarized data. In what follows, we shall use < a/ < 0.5, the range proposed in Ref. 
[50] in agreement with Ref. [48]. 

We have to generalize this scenario to describe off-forward phenomena. Of course, the 
forward limit of our GPDs formula, Eq. (38), has to be given by Eq. (57). By taking the 
forward limit of Eq. (38), one obtains: 

H q (x,0,0)=T, -H qo (z,0,0)H qoq (-,0,0) 

qo Jx Z \Z J 

= -^0(^(7,0,0) , (74) 

qo Jx Z \Z / 

so that, in order for the latter to coincide with Eq. (57), one must have H qoq (x, 0, 0) = 
A(f) qoq (x). 

In such a way, through the ACMP prescription, the forward limit of the unknown con- 
stituent quark GPD H qoq (^ , J, A 2 ) can be fixed. 

The off-forward behavior of the Constituent Quark GPDs has to be modelled and it will 
be done, as in Ref. [29] for the unpolarized case, using the "a-Double Distributions" (DD's) 
language proposed by Radyushkin [36,37]. The DD's, a, A 2 ), are a representation of 
GPDs which automatically guarantees the polynomiality property. GPDs can be obtained 
from DD's after a proper integration, and, as is explained in [36,37], parton densities are 
recovered in the forward limit, while meson distribution amplitudes are obtained in the P = 
limit of the DD's. In some cases, such a transparent physical interpretation, together with 
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the symmetry properties which are typical of distribution amplitudes (a — > —a symmetry), 
allows a direct modelling, already developed in [16]. 

The relation between any GPD H, defined a la Ji, for example the one we need, i.e. H qoq 
for the constituent quark target, is related to the a-DD's, which we call A$ qoq (x, a, A 2 ) for 
the constituent quark, in the following way [36,37]: 



fl fl— \x\ 

H qoq (x, f , A 2 ) = / dx 5{x + £a- x)A® qoq ( x, a, A 2 )da . (75) 

J-l J-l+\x\ 



With some care, the expression above can be integrated over x and the result is explicitly 
given in [37], where a factorized ansatz is suggested for the DD's which, for the helicity- 
dependent one, could be written: 

A$ qoq (x, a, A 2 ) = h(x, a)A$, 0f ,(5)G«,(A 2 ) , (76) 

with the a dependent term, h(x,a), having the character of a mesonic amplitude. 

Besides, in Eq. (76) A<& qoq (x) represents the forward density and, eventually, G qo (A 2 ) 
the constituent quark axial ff. 

One immediately realizes that the helicity-dependent GPD of the constituent quark, Eq. 
(75), with the factorized form Eq. (76), fulfills the crucial constraints of GPDs, i.e., the for- 
ward limit, the first-moment and the polynomiality condition, the latter being automatically 
verified in the DD's description. 

In the following, it will be assumed for the constituent quark GPD the above factorized 
form, so that one needs to model the three functions appearing in Eq. (76), according to 
the present description. 

For the amplitude h, use will be made of one of the simple normalized forms suggested 
in [36], on the bases of the symmetry properties of DD's: 

^ 3(l-5) 2 -a 2 
= 4 (l-> • ^ 

In addition, since we will identify quarks for x > £/2, pairs for x < |£/2|, antiquarks for 
x < — £/2, and, since in our approach the forward densities AQ qoq (x) have to be given by 
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the standard A$ functions of the AC MP approach, Eqs. (71)-(73), one has, for the DD of 

flavor q of the constituent quark: 

_ r(/ i (5,a)i» (? ^(5) + / i (5,a)<l^(x))G;(A 2 ) for x > 

A% )q (x,a,A 2 ) = { _ (78) 

I -^(-^,«)^ 093 (-x)^ o (A 2 ) for 5 < 

The above definition, when integrated over a gives the correct limits [37]: 

f 1 X A$ qoq (x, a, A 2 = 0)U >0 da = A<f> goq (x) , (79) 

J-l+X 

and 



f 1 A$ qoq (x, a, A 2 = 0)U <0 da = -A$ qoq (-x) . (80) 

J-l+\x\ 



Eventually, as axial ff we will take a simple monopole form corresponding to a constituent 
quark size tq ~ 0.3/m: 

G qo (A 2 ) = , (81) 

showing i.e. the same Q 2 dependence of the electromagnetic ff used in [29], the latter 
scenario being strongly supported by the analysis of [33]. 

By using such a ff and Eq. (77), together with the standard ACMP A$'s, Eqs. (71) - 
(73), in Eq. (78), and inserting the obtained A§ qoq (x, a, A 2 ) into Eq. (75), the constituent 
quark GPD in the ACMP scenario can be eventually calculated. 

Let us discuss our factorized choice for the DD of the constituent quark, Eq. (76). Such 
a factorization does not hold for the GPD of the proton [4]. No microscopic calculation 
supports it [23]. The present approach is consistent with the absence of factorization of the 
proton GPD in A 2 — dependent and A 2 — independent terms since an additional, crucial, A 2 
dependence comes from the quantity H qo in Eq. (38), i.e., from the proton wave function of 
the quark model used. 

The factorized form is assumed for the GPD of the constituent quark for simplicity. We 
do not know what the helicity-dependent GPD of the constituent quark is and we start 
modelling it by using the simplest possible structure which fulfills the required constraints. 
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Clearly this choice deserves a deeper explanation, which we do not have at the moment. 
However, one expects that the structure of the constituent quark GPD is closer to the 
asymptotic behavior, than that of the proton, because it probes higher energies. The non 
perturbative contributions arise in the precise shape of the three functions building the 
GPD, which we take from general QCD based arguments and from data. It is true that even 
the forward structure of the constituent quark is not known, but we use a parameterization 
which performs well in describing the data [31]. In this sense, the forward limit is well 
reproduced, as shown in the next section. On the contrary, the profile Eq. (77) is assumed 
to be valid also for the constituent quark, even in the helicity dependent case, simply because 
it satisfies some of the required properties. 

These arguments and a small set of parameters, fixed once forever, make our approach 
able, if a proper CQM is used, to describe DIS data for parton distributions, elastic form 
factors, and to predict unpolarized and polarized GPDs. It allows for an effective description 
of the known structure of the nucleon and its predictive power needs to be tested with new 
observables. 

V. RESULTS AND DISCUSSION 

In this section we present the results obtained for the helicity-dependent GPD H q of the 
proton, for ^ 2 C 1 and A 2 <C m 2 , according to the approach described so far. The main 
equation to be evaluated is Eq. (38), written again here below for the sake of clarity: 

qo Jx Z \Z Z ) 

In the above equation, the quantity H qoq , the helicity-dependent constituent quark GPD, 
is modelled according to the arguments described in the previous section. This means that 
it is obtained evaluating Eq. (75), where the DD of the constituent quark, A<& qoq (x, a, A 2 ) 
, is given by Eq. (78), calculated in turn through the ff, Eq. (81), the function h, Eq. (77), 
together with the forward ACMP A$'s, Eqs. (71) - (73). 
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The other ingredient in Eq. (38), H qo , has been evaluated according to Eq. (22): 

H qo (z, £, A 2 ) = J dpn qo (p, p + A)5^z + ^-t- 

The calculation has been performed in the Breit frame, where one has, in the NR limit 
studied, \/2P + — > M. The off-diagonal momentum distribution appearing in the formula 
above, n qo (p, p + A), defined in Eq. (44), has been evaluated within the Isgur and Karl (IK) 
model [24] , according to the definition 

n qo (p,P + A) = n qo+ (p,p + A) - n qo _(p,p + A) , (82) 

— # 

where n qo± (p,p + A) is the off-diagonal momentum distribution for the valence quark 
momentum p and spin projection parallel or antiparallel to that of the nucleon. n q ±(p, p+ A) 
can be evaluated projecting out the appropriate spin and flavor component of the constituent 
quark and in the corresponding baryonic state is given by 

n Uo /do,±ip,P + A) = 3 ( s/ < N, J z = +I]i^p_i^p|7V, J z = +1 > s/ ) x 

x J dp 1 dp 2 ip*(pi,P2,p + A)^(pi,P2,p) , (83) 

where \N, J z = +| > s / is the spin-flavor state and ip{pi,P2,Ps) is the proton wave function 
in momentum space (cf. Ref. [23] for the definition in the unpolarized case). 

In the IK CQM [24], including contributions up to the 2hu shell, the proton state is 
given by the following admixture of states 

\N) = a s \ 2 S 1/2 )s + a S '\ 2 S[ /2 )s + a M \ 2 S 1/2 )M + a v \ 4 D 1/2 ) M , (84) 

where the spectroscopic notation \ 2S+1 Xj) t , with t = A,M,S being the symmetry type, 
has been used. The coefficients were determined by spectroscopic properties to be [51]: 
a s = 0.931, a s , = -0.274, a M = -0.233, a v = -0.067. 

The results for the GPD H(x, £, A 2 ), neglecting in (84) the small D-w&ve contribution, are 
found to be: 




27 



M ( 3 \ 3 ^ 2 a 2 r r 
H u (x, f , A 2 ) = 3— ( — ) / dk x / dk y f (k x , k y , x, £, A 2 ) 



x 



a 3 \2%J 

Af s (k x , k y , x, £, A 2 ) + Af(k x , k y , x, £, A 2 ) 



(85) 



iW / 3 \ 3 ^ 2 a 2 f f 
H d (x, £, A 2 ) = 3— ( — ) J dk x J dk y f (k x , k y , x, f , A 



x 



a 3 \2nJ 

-Af s (k x , k y , x, £, A 2 ) + \Af{k x , k y , x, i, A 2 ) + f(k x , k y , x, f , A 2 ) 



for the flavors it and d, respectively, with 



fo(k x ,k y ,x,£,A ) 



fc + k 2 



fa{A X , k X )f a (Ay, ky)f a (A z , k Z ) , 



, (86) 



(87) 



/ a (A^)=e-^U 



(88) 



M 2 (a; + 2 -(m 2 + fcg + fcg) 
2M(x + 



(89) 



4 4 

A/ a (fc x , fey, x, £, A 2 ) = -a 2 + -a 2 



5 3 A; 2 3 /c 4 1 / A 2 ^ A /A; 2 _ N 
4 2a 2 + 4a 4 + a 2 g ~*~ j L 2 . 



?-^ + ^JLf_4! + A.^-i N 



4 a 2 8 a 4 2a 2 \ 3 



-1 + 



2 /A 5 



3a 2 V 3 



- A - k 



(90) 



Af(k x ,k y ,x,t,A 2 ) = - 



+ 



1 5 
^23 

2e_ ( 

3a 4 V 



a M as> 



\_1f_ _2_ 

Q~a 2 + 2a* 3a 2 



/ A 2 - ^ 
A., 



A 2 - 



(91) 



/(A*, fcj,, x, £, A 2 ) = l - ^\k 2 - A 2 + 3A ■ k , (92) 

and k = ^m 2 + k 2 + k 2 + k 2 , being m ~ M/3 the constituent quark mass. Here we have 
used the notation k 2 = k 2 ,k = (k x ,k y ,k z ). Evaluating the forward limit, A 2 = 0, £ = 0, 
of Eqs. (85) and (86), the expressions of the polarized quark densities in the IK model, 
obtained in [41], are consistently recovered. 
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The harmonic oscillator parameter, a, of the IK model, has ben chosen so that the 
experimental r.m.s. of the proton is reproduced by the slope, at A 2 = 0, of the charge ff, so 
that a value of a — 1.18/ -1 is obtained [29]. Such a ff reproduces well the data at the low 
values of A 2 which are accessible in the present approach. For higher values of A 2 , it would 
not be realistic [51]. The same quality of results is expected for the proton axial form factor, 
which is obtained by integrating Eq. (38) over x, according to Eq. (42). The experimental 
dipole trend [52] is qualitatively recovered at the low values of A 2 under scrutiny here. 

Results for the w-quark H u distribution, at the scale of the model /i 2 , are shown in Figs 

2 to 4. In Fig. 2, it is shown for A 2 = —0.1 GeV 2 and £ = 0.1. One should remember that 
the present approach does not allow to estimate realistically the region —A 2 > m 2 ~ 0.1 
GeV 2 , so that we are showing here the result corresponding to the highest possible A 2 
value. Accordingly, the maximum value of the skewedness is therefore £ ~ 0.17 (cf. Eq. 
(2)), fulfilling the requirement £ 2 <C 1. The dashed curve represents what is obtained in the 
pure IK model, i.e., by evaluating Eq. (85). One should notice that such a result, obtained 
in a pure valence CQM, should vanish for x < £ and for x > 1. The small tails which are 
found in these forbidden regions represent the amount of support violation of the approach. 
In particular, for the shown values of A 2 and £, a violation of 2 % is found. In general, in 
the accessible region the violation is never larger than few percents. The full curves in Fig. 

3 represents the complete result of the present approach, i.e., the evaluation of Eq. (38) 
following the steps and using the ingredients described in this section and in the previous 
one. The two curves are obtained in correspondence of the two extreme values of the Regge 
intercept [48], a/ = (upper curve) and a/ = 0.5 (lower curve), used to define the helicity 
dependent constituent quark structure functions [31]. A relevant contribution is found to 
lie in the ERBL region, in agreement with other estimates [10]. As already explained, the 
knowledge of GPDs in the ERBL region is a crucial prerequisite for the calculation of all 
the cross-sections and the observables measured in the processes where GPDs contribute. 
We notice that the ERBL region is accessed here, with respect to the approach of Ref. 
[23] generalized to the helicity-dependent case, which gives the dashed curve, thanks to the 
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constituent structure which has been introduced implementing the ACMP procedure. A 
similar result has been found in [29] in the unpolarized case. 

In Fig. 3, special emphasis is devoted to show the ^-dependence of the results. For the 
allowed £ values, H u (x, £, A 2 ), evaluated using our main formula, Eq. (38), is shown for 
four different values of x. It is clearly seen that such a dependence is stronger in the ERBL 
region, than in the DGLAP region, in good agreement with other estimates [10] and with 
our previous analysis, performed in the unpolarized case [29]. To allow for a complete view 
of the outcome of our approach, in Fig. 4 the x and £ dependences are shown together in a 
3-dimensional plot. 

The results shown so far are associated with the low scale of the model, the hadronic 
scale nl, fixed to the value 0.34 GeV 2 , as discussed in Section 4. 

In order to have a feeling of the quality of the present predictions, since no data are 
available at the moment for H(x, £, A 2 ), we evaluate in the following the forward limit 
of the obtained quantities. In Figs. 5 and 6, the proton and neutron spin dependent 
structure function gi(x,Q 2 ), built from the forward limit of Eq. (38) evaluated for q = 
q v ,q s ,9, are compared with experimental data. In Ref. [31], the succesful analysis had been 
performed using the model of Ref. [53] , so that the encouraging results shown there have to be 
confirmed when the IK model is considered, as it is being done here. To evaluate g\(x,Q 2 ) 
at the experimental scale, a NLO evolution of the model parton distributions has been 
performed. It is known that perturbative QCD to this order allows the proposal of varied 
factorization schemes [54] , whereby one is able to define the partonic distributions in different 
ways without altering the physical observables. In our analysis the parton distributions are 
determined by the quark model through the ACMP procedure. It is evident that different 
factorization schemes lead to different results. Following Ref. [31], we have adopted on 
physical grounds the AB scheme as defined in ref. [55] 1 . In it some relevant physical 



It consists in modifying minimally the NLO polarized anomalous dimensions [56], calculated in 
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observables, such as AE, the total spin carried by the quarks, are scale independent, i.e., they 
behave as conserved quantities, and therefore the partons have a well defined interpretation 
in terms of constituents. Fig. 5 and 6 show that the present approach describes well the 
data [60,61] for the proton structure function g±, and the trend of the neutron data [62] is 
also reasonably reproduced. The quality of the agreement is comparable with that obtained 
in Ref. [31] using the model of [53]. 

The present approach in the forward limit, i.e. the description given in Ref. [31], also 
allows for the predictions of axial charges of the proton. By using the IK model, the following 
values are obtained (the first listed value corresponds to the Regge slope a/ = 0.5, the 
second to aj = 0.): for the octet charge, a 8 = 0.540 0.671 (experimental value: 0.58 ±0.03 
[58]), for the isospin charge, a 3 = 0.857 -h 1.063 (experimental value: 1.257 ± 0.003 [59]), 
while the predictions for the singlet charge and for the spin carried by the gluons lead to 
AS = 0.402 + 0.530 (from the analysis of Ref. [50] : 0.45 ± 0.09). Again, the agreement 
with data of these results, obtained by using the IK wave functions, is comparable to the 
one obtained in Ref. [31] by using the wave functions of Ref. [53]. 

This finding deserves a further comment. The lack of high momentum components in 
the IK wave function with respect to that of the model of Ref. [53] would indicate the latter 
as a better tool to describe the data. Nevertheless, the low-x behavior of the model results 
turns out to be mainly governed by the constituent quark structure rather than by the 
constituent quark wave function. The correct Regge-behavior given by our parametrization 
of the constituent quark structure functions provides us with a model independent prediction. 
Such a feature has not to do with the use of a relativistic model, which would improve only 
the high x behavior. As discussed in Ref [63], relativity and constituent quark structure are 
independent features, both of them necessary if model calculations are used for a proper 



the MS scheme, in order to have the axial anomaly governing the first moment of gi, as proposed 
in Ref [57]. 
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description of data in the full x range. 

Since no data are presently available for helicity-dependent GPDs, it is certainly useful to 
compare the results of model calculations with model independent constraints, if available. It 
has been proposed that they could be obtained from analysis of data in particular kinematical 
conditions. Very recently, this issue has been addressed in two papers, where GPDs have 
been obtained, at £ = 0, by fitting the first moment of a reasonable ansatz of GPDs to ff 
data [64,65]. Anyway, such a procedure is difficult to be realized for the quantity under 
scrutiny here, due to the lack of a satisfactory experimental knowledge of the axial ff, which 
cannot provide us with a clear constraint on H q (x, 0, A 2 ). As a matter of fact, in Ref. [65] no 
prediction has been produced for it, so that no comparison is possible between the present 
results and those of [65]. On the contrary, some curves reproducing H q (x, 0, A 2 ) have been 
shown in [64]. Nevertheless, as the authors of that paper clearly state, the lack of precise data 
for the axial ff does not allow a realistic parameterization of H q (x, 0, A 2 ), but rather a simple 
ansatz for it. This ansatz is given by the product of the polarized parton distributions, taken 
from the "scenario 1" of Ref. [66] at Q 2 = 4 GeV 2 , and an exponential term, giving the A 2 
dependence and additional x dependence, taking into account the different Regge behavior 
at different values of the momentum transfer. The latter dependence is borrowed from the 
analysis performed for the helicity independent GPDs, giving the so called default fit of [64], 
and no argument is given to support this procedure in the helicity dependent case. Anyway, 
we show in Fig. 7 the comparison of our results for the valence distribution H"(x,0, A 2 ), 
evaluated at A 2 = —0.1 GeV 2 , evolved to NLO from the hadronic scale up to Q 2 = 4 GeV 2 , 
with the ansatz proposed in [64], evaluated in the same kinematical conditions. It is seen 
that the latter and the present model predict a similar scenario. The very low x behavior is 
governed in both cases by the Regge phenomenology, which is left here free to vary between 
the two full lines of the figure, while it is fixed in Ref. [64] in the way described above, 
producing a big enhancement at low values of —A 2 . 

Summarizing this section, we have developed here a scheme which provides us with the 
helicity-dependent GPD H in the full x range. This is obtained thanks to the constituent 
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quark structure, implemented dressing the three quarks of a CQM, where initially only the 
DGLAP region of GPDs was accessible. This is an important development, a prerequisite for 
any attempt to calculate cross sections and asymmetries of related processes. The next step 
of our studies will be indeed to use the obtained GPDs for the evaluation of cross sections. 



VI. CONCLUSIONS 

A thorough analysis of both polarized and unpolarized data have shown, in previous 
work, that constituent quarks cannot be considered elementary when studied with high- 
energy probes. Certainly evolution is needed, i.e. the constituent quarks at the hadronic 
scale have to be undressed by incorporating bremsstrahlung in order to reach the Bjorken 
regime. Besides, the constituent quarks should be endowed with soft structure in order 
to approach the data. Thus, the constituent quarks, when under scrutiny by high-energy 
probes, seem to behave as complex systems, with a very different behavior from the current 
quarks of the basic theory. These features, which we found in structure functions, have also 
been recently discussed in form factors [33]. 

The formalism of composite constituent quarks has been recently applied to the study 
of generalized parton distributions, developing a formalism which expresses the hadronic 
GPDs in terms of constituent quarks GPDs by means of appropriate convolutions [29]. The 
approach has been extended here to treat helicity-dependent GPDs. To our knowledge, this 
represents the first calculation which gives, in a constituent quark framework, estimates 
of helicity-dependent GPDs in the ERBL region. Following the spirit of Ref. [29], this is 
done incorporating phenomenological features of various kinematical regimes. In particular, 
use has been made of the Radyushkin's factorization ansatz, thus our constituent quark 
GPDs are defined in terms of the product of three functions: i) the constituent quark 
polarized structure function, where we use the parameterization of [31] following the ACMP 
proposal [32]; ii) Radyushkin's double distributions [36]; iii) constituent quark form factor 
as suggested in Ref. [33]. Once the helicity-dependent GPDs are defined in this way, we 
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have developed the scheme to incorporate them into any nucleon model by appropriate 
convolution. In order to show the type of predictions to which our proposal leads, we have 
used here, as an illustration, the IK [24] model. However, in this latest step of our scheme, 
any non-relativistic (or relativized) model can be used to define the hadronic GPDs. 

As in the unpolarized case [29] , the present scheme transforms a hadronic model, in whose 
original description only valence quarks appear, into one containing all kinds of partons (i.e., 
quarks, antiquarks and gluons). Moreover, the starting model produces no structure in the 
ERBL region, while after the structure of the constituent quark has been incorporated, it 
does. The completeness of the x-range, for the allowed A 2 and £ ( A 2 << m 2 and ^ 2 <C 1 
), of the present description, is a prerequisite for the calculation of cross-sections and other 
observables in a wide kinematical range. Relativistic corrections, which permit to access a 
wider kinematical region, could be included in the approach. 

We have reminded the reader of the calculation for the forward polarized structure func- 
tions to see how in this case, where experimental data are available, our scheme leads, even 
with a naive quark model, to a reasonable description of the data. Thereafter, we have 
proceeded to calculate the GPDs of physical interest to guide the preparation and analy- 
sis of future experiments. Our calculation seems to be consistent with an ansatz recently 
suggested in Ref. [64]. 

This work is the continuation of an attempt to describe the properties of hadrons in 
different kinematical and dynamical scenarios. Although our description can never be a 
substitute of Quantum Chromodynamics, it may serve to guide experimenters to physical 
processes where the theory might show interesting features, worthy of a more fundamental 
effort. 
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FIGURE CAPTIONS 



Fig. 1: The handbag contribution to the DVCS process in the present approach. 

Fig. 2: The GPD H for the flavor u, for A 2 = -0.1 GeV 2 and f = 0.1, at the momentum 
scale of the model. Dashed curve: result in the pure Isgur and Karl model, Eq. (85). The 
small tails which are found in the forbidden regions, x < £ and x > 1, represent the amount 
of support violation of the approach. Full curve: the complete result of the present approach, 
Eq. (38). The two curves are obtained in correspondence of the two extreme values of the 
Regge intercept [48], a/ = (upper curve at low x) and a/ = 0.5. 

Fig. 3: For the £ values which are allowed at A 2 = —0.1 GeV 2 , H u (x, £, A 2 ), evaluated using 
our main equation, Eq. (38), is shown for four different values of x, at the momentum scale 
of the model. From top to bottom, the dash-dotted line represents the GPD at x — 0.05, 
the full line at x = 0.1, the dashed line at x — 0.2, and the long-dashed line at x — 0.4. 
The curves are obtained in correspondence of the value of the Regge intercept a/ = 0.5 (see 
text). 

Fig. 4: The x and £ dependences of H u (x, £, A 2 ), for A 2 = —0.1 GeV 2 , at the momentum 
scale of the model. The value of the Regge intercept has been taken to be aj = 0.5 (see 
text). 

Fig. 5: We show the structure function xg[(x, Q 2 ) obtained at Q 2 = 10 GeV 2 by evolving 
at NLO the model calculation [24] considering the structure of the constituents, for the two 
extreme Regge behaviors mentioned in the text (a/ = is the upper curve, here and also in 
the following figure). The data from refs. [60] (full dots) and [61] at Q 2 10 GeV 2 are also 
shown. 

Fig. 6: The structure function xg™{x, Q 2 ) for the neutron evolved at NLO to Q 2 = 5 GeV 2 , 
for the two extreme Regge behaviors mentioned in the text, is shown by the two full curves. 
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The data from refs. [62] at Q 2 = 5 GeV 2 are also shown. 

Fig. 7: Comparison of the results of the present approach for the valence distribution 
H^{x, 0, A 2 ) (full curves), evaluated at A 2 = —0.1 GeV 2 , evolved to NLO from the hadronic 
scale up to Q 2 = 4 GeV 2 , with the ansatz proposed in [64] (dotted curve), evaluated in the 
same kinematical conditions. The two full curves are obtained in correspondence of the two 
extreme values of the Regge intercept [48], a/ = (upper curve at low x) and a/ = 0.5. 
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